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We consider the continuous-time version of our recently proposed quantum theory of optical 
temporal phase and instantaneous frequency [Tsang, Shapiro, and Lloyd, Phys. Rev. A 78, 053820 
(2008)]. Using a state- variable approach to estimation, we design homodyne phase-locked loops that 
can measure the temporal phase with quantum-limited accuracy. We show that post-processing can 
further improve the estimation performance, if delay is allowed in the estimation. We also investigate 
the fundamental uncertainties in the simultaneous estimation of harmonic-oscillator position and 
momentum via continuous optical phase measurements from the classical estimation theory perspec- 
tive. In the case of delayed estimation, we find that the inferred uncertainty product can drop below 
that allowed by the Heisenberg uncertainty relation. Although this result seems counter-intuitive, 
we axgue that it does not violate any basic principle of quantum mechanics. 
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I. INTRODUCTION 

Optical phase measurements at the fundamental quan- 
tum limit of accuracy are an important goal in science 
and engineering and crucial for future metrology, sens- 
ing, and communication applications. While the single- 
mode case has been extensively studied, less attention 
has been given to the quantum measurements of a tem- 
porally varying phase. Theoretically, the temporal-phase 
positive operator-valued measure (POVM) describes the 
optimal quantum measurements pj], but it is difficult 
to perform such measurements in practice. Adaptive 
homodyne detection [3, Q is a much more feasible ap- 
proach, and Berry and Wiseman have proposed the use 
of a homodyne phase-locked loop to estimate the phase 
when the mean phase is a classical Wiener random pro- 
cess d, |5[. On the other hand, we have recently shown 
in Ref. [l|] how homodyne phase-locked loops can be 
designed using classical estimation theory to perform 
quantum-limited temporal phase measurements when the 
mean phase is any stationary Gaussian random process. 

The main purpose of this paper is to unify and general- 
ize the two distinct approaches undertaken by Berry and 
Wiseman and ourselves, under the common framework 
of classical estimation theory. In Sec. |TT1 we first extend 
our discrete-time theory proposed in Ref. jlj to the con- 
tinuous time domain. In Sec. IIII Al we generalize Berry 
and Wiseman's results to a much wider class of random 
processes using the Kalman-Bucy filtering theory [gl, 0] ■ 
The Kalman-Bucy approach guarantees the real-time es- 
timation efficiency provided that the phase-locked loop 
operates in the linear regime. Our approach also signif- 
icantly simplifies the design of phase-locked loops, com- 
pared to the more computationally expensive Bayesian 
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state estimation approach suggested by Pope et al. [1]. 
In Sec. IIII B1 we show that the Wiener filtering technique 
used in our previous paper [l[ is equivalent to Kalman- 
Bucy filtering at steady state. In Sec. IIVI we point out 
that Berry and Wiseman's results are not optimal if delay 
is permitted in the phase estimation process, and post- 
processing can further improve the phase estimation per- 
formance beyond that offered by Kalman-Bucy or Wiener 
filtering. We illustrate these concepts by considering 
the specific cases of the mean phase being an Ornstein- 
Uhlcnbeck random process as well as the Wiener process 
studied by Berry and Wiseman. Apart from the theoret- 
ical importance of our results in the context of quantum 
estimation and control theory, they should also be of im- 
mediate interest to experimentalists and engineers who 
wish to achieve quantum-limited temporal phase mea- 
surements, as we expect our proposals to be realizable 
using current technology. 



In Sec. |Vl we investigate the fundamental problem of 
simultaneous harmonic-oscillator position and momen- 
tum estimation at the quantum limit by continuous op- 
tical phase measurements. The problem can be cast di- 
rectly in the framework of classical estimation theory for 
Gaussian states. The use of Kalman-Bucy filtering for 
real-time position and momentum estimation has been 
proposed by Belavkin and Staszewski [8] and Doherty et 
al. [9| , who have shown that the quantum state of the har- 
monic oscillator conditioned upon the real-time measure- 
ment record is a pure Gaussian state. Here we show that 
the inferred position and momentum estimation errors 
according to classical estimation theory can be further 
reduced below the Heisenberg uncertainty product, if de- 
lay is allowed in the estimation. While counter-intuitive, 
we explain in Sec. IV CI why this result does not violate 
the basic principles of quantum mechanics. 
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II. PHASE IN THE CONTINUOUS TIME 
DOMAIN 

For completeness, we first review the continuous time 
limit of our discrete-time theory of temporal phase [H, 
as previously described in Ref. ■ Consider the optical 
envelope annihilation and creation operators A{t) and 
A^{t), respectively, in the slowly varying envelope regime, 
with the time-domain commutation relation 



i(i),it(t') =6{t-t') 



(2.1) 



Let dn{t) be a continuous- time discrete-photon- number 
random process, a realization of which is depicted in 
Fig. [1] and tj be the times at which dn{Tj) is non-zero. 
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FIG. 1; (Color online). A realization of the continuous-time 
discrete-photon-number random process. 

A Fock state with a definite dn{t) can be defined as 



\dnit)) 



n 



A\Tj)Vdi 



dn{Tj ) 



|0), (2.2) 



which is an eigenstate of the photon-number flux opera- 
tor A^{t)A{t), 



A^{t)A{t)\dn{t)) ^I{t)\dn{t)), (2.3) 
=Y,dn{T,)S{t-r,). (2.4) 



The Fock states form a complete orthogonal basis of the 
continuous- time Hilbert space. 



J2 \dnit)){dnit)\=l 

dn(t) 



(2.5) 



where the sum is over all realizations of dn(t). For a 
quantum state /5, the photon-number probability distri- 
bution is 

P[dn{t)] = TY{p\dn{t)){dn{t)\} , ^ P[dn{t)] = 1. 

dn{t) 

(2.6) 



For example, a coherent state is defined as 



N 



dtA{t)A^{t) 



|0), 



\A{t)) = exp 

N=l dt\A{t)\\ A{t)\A{t)) =A{t)\A{t)), 



(2.7) 



where A{t) is the mean field. The photon-number prob- 
ability density is then 



P[dn{t)\ = lim e 

St^dt 



-N 



n 



[\A{t,)\Ht_ 



dn{t,)\ 



tj = to + jSt, 



(2., 



which describes a Poisson process, as is well known 

A temporal phase state can be defined as the functional 
Fourier transform of the Fock states. 



dn(t) 



\^dn{Tj)4){Tj 



r poo 
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\dn{t)) 

\dn{t)). (2.9) 



In terms of the temporal phase states, a temporal-phase 
POVM can be defined as 



\im)]^\mm)i 



(2.10) 



which is the continuous limit of the one defined in Ref. [l| 
and can be normalized using a path integral with the 
paths restricted to a range of 27r, 



/ Dcl>[t)tl m)] - i, Dc^it) EE lim TT 



2t: ' 



0o(i)<0(t)<0o(O + 2^- (2.11) 
The temporal-phase probability density is thus given by 

pm)] = Tr [mm]} , J mmm] = i- (2.12) 

It is difficult to analytically calculate p[(/)(t)] for most 
quantum states of interest, so perturbative or numerical 
methods should be sought. 

For the design of homodyne phase-locked loops, the 
Wigner distribution is of more interest. For a Gaussian 
state with uncorrelated quadratures, it can be written as 

\Y. I dtdT^,{t)K-\t,T)i,{T) , (2.13) 



(X exp 
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where (t) are quadrature processes, 

- (^A{t)e-'^'^'^ + A*(^)e'''(*)^ , (2.14) 

- {~i [A(t)e~*''(*) - yl*(t)e*''(*)] ^ , (2.15) 

A{t) is the complex field variable in phase space, 9 is an 
arbitrary phase, and KJ^it^r) is defined in terms of the 
covariance functions Kjit^r) as 

j duKj{t,u)K-^{u,T) = S{t - t), (2.16) 
K,{t,r)^{^,m,{T)). (2.17) 

The covariance functions must satisfy the uncertainty re- 
lation 



duKi{t,u)K2{u,T) > 6{t - r), 



(2.18) 



which becomes an equality if and only if the state is 
pure. In particular, the covariance functions for a co- 
herent state are 



Kl{t,T) ^ K2{t,T) =dit-T). 



(2.19) 



III. PHASE-LOCKED LOOP DESIGN 



A. Kalman-Bucy filtering 



a X 




I exp I 



FIG. 2: (Color online). A homodyne phase-locked loop. LO 
denotes local oscillator. 

Consider the homodyne phase-locked loop illustrated 
in Fig. [2] The output of the homodyne detection can be 
written as 



(3.1) 



where (f>{t) is the mean phase of the optical field, which 
contains the message to be estimated, 4>'{t) is the local- 
oscillator phase, and z{t) is the quantum noise. For a 



phase-squeezed state with squeezed quadrature ^2 (t) and 
anti-squeezed quadrature ^i{t), z{t) can be written as 

+ Ut)cosm)-(t>'{t)]}, (3.2) 

where A = (A) = |^| exp(z0) is the mean field. For 
generality, we let the message be a vector of n random 
processes, 



xit) 



Xl{t) 



(3.3) 



with the mean phase proportional to the first one, 

^{t) = Cit)x{t), C(t) = [/3,0,...,0]. (3.4) 

In the Kalman-Bucy formalism, x{t) is modeled as zero- 
mean random processes that satisfy a set of linear differ- 
ential equations. 



dx(t) 
dt 



A{t)x{t) +B{t)u{t), 



(3.5) 



where A{t) and B{t) are n x n and n x m matrices, 
respectively, and u(t) is a vector of m zero-mean white 
Gaussian inputs with autocorrelation 



{u{t) (g) M(r)) US{t - t). 



(3.6) 



We focus on coherent states, so that z{t) can be modeled 
as an independent white Gaussian noise according to its 
Wigner distribution. 



(z(i)z(T)> -Z(t)J(i-r), Z(t) 



4|^|^ 



4^' 



(3.7) 



where ujq is the optical carrier frequency and V is the av- 
erage optical power. The additive white Gaussian noise 
allows us to apply classical estimation techniques directly. 
Coherent states should also be of more immediate inter- 
est to experimentalists and engineers, as they are easier 
to generate and more robust to loss compared to non- 
classical states. For a phase-squeezed state, the statistics 
of z{t) depend on 0(i) — (p'{t), but one may still wish to 
approximate z(t) as an independent Gaussian noise by 
neglecting the anti-squeezed quadrature ^i{t), in order 
to take advantage of classical estimation techniques 

The purpose of the phase-locked loop is to make (j)' (t) 
the optimal estimate of 4>{t), using the measurement 
record of ri{T) in the period to < t < t, such that we 
can linearize Eq. (|3.ip . 



7j{t)^^t)-q^'{t)+zit), 



(3.8) 



when the following condition, called the threshold con- 
straint in classical estimation theory P, @], is satisfied. 



([0(O-0'(O]')«1- 



(3.9) 
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The threshold constraint ensures that the phase-locked 
loop is phase-locked. 

If the canonical measurements characterized by the 
temporal-phase POVM can be performed, we can instead 
modulate the phase of the incoming field by —<j>'{t) and 
perform the canonical measurements, producing an out- 
put 



r,cit)^f{Ht)-^'it)+^(t)), 



(3.10) 



where /((^) must be a periodic function, such as a saw- 
tooth function, 



/((/)) = [((/) — tt) mod 27r] — tt. 



(3.11) 



and z(t) is the quantum phase noise and independent of 
(/)(t) and (t)'{t) for any quantum state. Because (/)(i) may 
exceed the 2tt range, it is still necessary to use the phase- 
locked loop to perform phase unwrapping. The follow- 
ing analysis can be applied to canonical temporal-phase 
measurements and arbitrary quantum states if ric{t) is 
linearized as 



(3.12) 



and z(t) is approximated as a white Gaussian noise. The 
same threshold constraint given by Eq. (|3.9p ensures that 
the periodic nature of r]c(t) can be neglected and the 
linearization is valid. 

The linearization allows us to use Kalman-Bucy filter- 
ing to produce the real-time minimum-mean-square-error 
estimates of x{t) ^B,^, which we denote as x'{t), 



dx' , 
—— = Ax 

dt 



rr?. 



(3.13) 



This is called the Kalman-Bucy estimator equation, r] 
is called the innovation, defined in terms of a general 
vectoral observation process y{t) as 



yit) = C(t)x{t) + z{t), r]{t) 



y{t) - C{t)x'{t), 

(3.14) 



where z{t) is a vectoral Gaussian white noise with mean 
{z{t)) = and covariance {z{t) » z{t)) = Z{t)5{t - r). 
For phase-locked loops, the homodyne output rj{t) can 
be used directly as the innovation, so z{t) = z{t) and 
Z{t) = Z{t). r is called the gain, given by 



hujo 



S21 



(3.15) 



and S is the estimation covariance matrix, defined as 

S(t) = {[x{t) - x'it)] ® [x{t) - x'it)]) , (3.16) 
which satisfies the variance equation, 

— = A'E + 'SA^ -■SC'^Z-^C'S + BUB'^. (3.17) 
dt 



Equations (|3.13p to (I3.17P are much simpler to solve than 
the conditional probability density equation suggested by 
Pope et al. for phase estimation The threshold con- 
straint becomes 



« 1, 



(3.18) 



and the initial conditions arc 



x'{t^) - {xit^)) = 0, S(io) = {x{to) ® x{to)). (3.19) 

Apart from phase estimation, Kalman-Bucy filtering can 
also be used to simultaneously estimate other parame- 
ters that depend linearly on the phase. The instanta- 
neous frequency, for instance, can be estimated by defin- 
ing X2 oc dxi/dt. The phase-locked-loop implementation 
of Kalman-Bucy filtering for general angle demodulation 
is depicted in Fig. [31 
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FIG. 3; (Color online). A phase-locked loop that implements 
Kalman-Bucy filtering for angle demodulation. 

For example, consider the message as an Ornstein- 
Uhlenbeck process. 



— — —kx -f Bu. 
dt 

The variance equation becomes 

dt hjjQ 
and the gain is 



(3.20) 



K = B^U, (3.21) 



T{t) = ^S(t). 
The variance equation can be solved analytically, 
M-^exph27(i-to)] 



^ -I- exp[-27(i - to)] 



_ 7/fc + 1 + AS(to) ^./^^i. 
- 7/fc-l-AS(to)' ^~ \k ^ 



A = 



hwok 



(3.22) 

(3.23) 

1/2 

(3.24) 



where the subscript ss denotes the steady state. 



1 

A 



kA 



1/2 



t - to > 



1 



/ K N kA 

[kA/ ' ~ 



» 1, 



7 

(3.25) 
(3.26) 
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and the threshold constraint is 

When the message is a Wiener random process, 
dx 



dt 



Bu, 



(3.27) 



(3.28) 



we can either foUow the same procedure as before to de- 
rive the Kalman-Bucy filter, or take the results for the 
Ornstein-Uhlenbeck process to the limit fc — > 0. Either 
way, assuming /3 = 1 for simplicity, we find 

^ + cxp[— 27(t - to)\ nujo 

(3.29) 



(3.30) 



^ _ 1 + 2ViVI](^o) 



7 = 2kVN, N = 



fkUnK 



At steady state, 

2VN 7 

(3.31) 



The threshold constraint is 



AN = > 1. 

TVjJqK 



(3.32) 



These results for the Wiener process agree with Berry 
and Wiseman's 



B. Wiener filtering 

In addition to the Kalman-Bucy state-variable ap- 
proach, Wiener's frequency-domain approach can also be 
used to design the phase-locked loop [l|, 0, [l^] ■ Defining 



y{i) = ^{t) + z{t), 



(3.33) 



it can be shown that Kalman-Bucy filtering is equivalent 
to the integral equation 



x'{t)= f dTH{t,T)y{T), 



(3.34) 



where H(t, r) is called the optimum realizable filter and 
satisfies the integral equation 

K^y{t,<j)= f dTH{t,T)Ky{T,a), (3.35) 

Jto 

K,y{t, a) = {x{t)y{a)) , Kyir, a) = {y{T)y{a)) . 

(3.36) 



If x{t) and y{t) are stationary and we let to — oo, 
Eq. (|3.35p becomes the Wiener-Hopf equation, 

K,yit-a)= j dTH{t~T)Ky{T-a), (3.37) 

which can be solved by a well-known frequency-domain 
technique [i|,[SIIl- 

For example, if x{t) is an Ornstein- 
Uhlenbeck process, its power spectral density in the limit 
of to ~oo is 



K,{t, a) = {x(t)x{a)) = K,{t - a), 
Sxi^^) — I dtKx{t) cxp{—iujt) = 

J-oo ^ - 

The power spectral density for y{t) is then 



fc2- 



Sy{uj) 



dtKy{t) exp(— iwi) 



To solve for H{t — t), we rewrite Sy{uj) as 



(3.38) 
(3.39) 



(3.40) 



\ 1/2 • , 

nwo \ luj -\- ^ 



AV 



iuj + k' 
(3.41) 



where 7 is given in Eq. (I3.24|) . and H+{lo) and \/H+{uj) 
are causal filters. Defining 



dtKxyit) exp(— iwt) 



cj2 + p 

the Wiener filter in the frequency domain is 



1 



, (3.42) 



(3.43) 



where the subscript -I- denotes the realizable part. To 
calculate the realizable part, first perform the inverse 
Fourier transform, 



doj Sxyjuj) 

2^i?;(w) 

1/2 



exp{iLL!t) 



■ktl 



(3.44) 



where U{t) is the Heaviside step function. The realiz- 
able part is then obtained by multiplying Eq. (13.441) by 
U{t) and performing the Fourier transform. After some 
algebra, we obtain 



H{co) 



icj -f 7 



7 — fc 



(3.45) 



To implement the Wiener filter in the phase-locked loop 
shown in Fig. 21 the loop filter that relates the homodyne 
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output 77 (t) to the estimate x'{t) is 



IV. SMOOTHING 



x'{t) = / dTL{t ~ T)r]{T) 



(3.46) 



dTL{t-r)[y{r)-f3x'{T)l (3.47) 



- H{uj), L{uj) - 



1 + l3L{iu) 



1 - l3H{uj) iuj + k' 
(3.48) 



The resuhing phase-locked-loop structure is equivalent to 
that obtained by Kalman-Bucy filtering at steady state, 
as both approaches implement the optimum realizable 
filter. 
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FIG. 4: (Color online). A phase-locked loop implementa- 
tion of Wiener filtering when the mean phase is an Ornstein- 
Uhlenbeck process. 



The mean-square error of Wiener filtering is given by 
the well-known expression [l|, 0, [l^ 



E 



hujQ 



duj 



In 



hujo 



1 

A 



kA 



1/2 
ll -1 



(3.49) 



which obviously must be the same as the steady-state 
error Sgs obtained by Kalman-Bucy filtering. The inter- 
ested reader is referred to Ref. ^ for an excellent treat- 
ment of Wiener filters. 

The advantage of Kalman-Bucy filtering over Wiener 
filtering is that the former can also deal with a wide class 
of nonstationary random processes that can be described 
by a system of linear equations (|3.5p , whereas Wiener fil- 
tering works only for stationary processes. In the special 
case of a Wiener process, however, we can first design a 
Wiener filter for an Ornstein-Uhlenbeck process and take 
the limit fc ^ 0. The result for /3 = 1 is 



r 

ss 



2kVN, 



1 



2VN 



(3.50) 



which is again the same as the steady-state Kalman-Bucy 
filter. 



Both Kalman-Bucy filtering and Wiener filtering pro- 
vide real-time estimates of x{t) based on the measure- 
ment record up to time t = t. If we allow delay in the 
estimation, we can use the additional information from 
more advanced measurements to improve upon the esti- 
mation. In the following we consider the optimal esti- 
mation of x{t) given the full measurement record in the 
interval to < t < T, also called smoothing in classical 
estimation theory 0, 0] . 



A. State-variable approach 

Given the output x' of the homodyne phase-locked 
loop designed by Kalman-Bucy filtering and the asso- 
ciated covariance matrix S, the optimal smoothing esti- 
mates of x{t), which we define as x{t), can be calculated 
using a state- variable approach, first suggested by Bryson 
and Frazier x{t) and the smoothing covariance 

matrix, 



n{t) EE {[x{t) - X{t)] (g> [X{t) - X{t)]) 



(4.1) 



can be obtained by solving the following equations back- 
ward in time. 



dx 



= Ax + BUB^H^^ {x - x') , 



-BUB^, 
with the final conditions, 
xiT) - x'iT), 



(4.2) 

T 



n + n{A + BUB^-E^ 

(4.3) 



n(T) = S(T). 



(4.4) 



In the io <C i <C T limit, we can calculate the steady-state 
smoothing covariance matrix Tigs by setting the right- 
hand side Eq. (|4.3p to zero and using the steady-state 
as S. 

Again using the Ornstein-Uhlenbeck process as an ex- 
ample, the steady-state smoothing error, also called the 
"irreducible" error , is given by 

(4.5) 



2fc(KA/fc-)- 1)1/2- 

This result is identical to that derived in [l[ using a 
frequency-domain approach. In the limit of A ;:!> fc/K, 



1 

2 VfcA/ 



(4.6) 



which is smaller than the error from Kalman-Bucy or 
Wiener filtering by approximately a factor of 2. For the 
Wiener process, the smoothing error is 



which is smaller than the filtering error by exactly a fac- 
tor of 2. 
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B. Two-filter smoothing 

An equivalent but more intuitive form of the optimal 
smoother was discovered by Mayne [l^ and Fraser and 
Potter [l^, who treat the smoother as a combination 
of two filters, one running forward in time to produce a 
prediction x' (t) via Kalman-Bucy filtering using the past 
measurement record, as specified by Eqs. p.l3|) to (|3.17|) . 
and one running backward in time to produce a retrodic- 
tion x"{t) using the advanced measurement record, 



dx" 

~dr 

'dt 



= Ax" Trj, 



= AS + SA^ 



with final conditions 



H-\T)a;"(T) = 0, H"'(r)=0. 



(4.8) 

BUB^, (4.9) 
(4.10) 



(4.11) 



The smoothing estimates and covariance matrix, taking 
into account both the prediction and the retrodiction, are 
given by 



5 = n (S-^a;' + H-^a;") 



n = (s 



-1 I TT-l"! 



(4.12) 
(4.13) 



The steady-state smoothing covariance matrix Tlss can 
be calculated by combining the steady-state predictive 
and retrodictive covariance matrices, 



C. Frequency-domain approach 



(4.14) 



For stationary Gaussian random processes and in the 
limit of to <^ t <^ T, a, frequency-domain approach can 
also be used to obtain the optimal smoother [l|, [l^] ■ 
The optimal smoothing estimates can be written in terms 
oiy{t) as 0,0, m 



xit) = [ dTG{t,T)y{T), 

Jtn 



where G{t, r) obeys 



K^y{t,c7)^ [ dTG{t,T)Ky{T,a) 

J to 



(4.15) 



(4.16) 



For K^y{t,a) = K^y{t - a), G{t,T) = G{t - r), and 
Ky{T, a) = Ky{T — (t), and in the limit of — ^ — oo and 
T ^ oo, we can solve Eq. (|4.16p by Fourier transform. 



G{lu) 



Sy{ui) ' 



(4.17) 



G{u}) is called the optimum unrealizable filter For an 
Ornstein-Uhlenbeck process, G{lu) is 



G(c.) 



4:f3V K 



(4.18) 



To implement this filter, one can use the homodyne 
phase-locked loop designed by Wiener filtering and a 
post-loop filter given by 



F(a;) 



G(w) fc + 7 



(4.19) 



(4.20) 



i/(ti>) —iio + 7 
The post-loop filter impulse response is 

/oo 7 
^F(c^)exp(*c.i) 

_ r (fc + 7)exp(7t), i < 0, 

~l 0, i>0, ^^■^^> 

which is anti-causal, so one must introduce a time delay 
td ^ 1/7 for F to be approximated by a causal filter. 
The optimal smoother designed by the frequency-domain 
approach is depicted in Fig. [5] 



% cxp i <j} ' 
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FIG. 5: (Color online). A homodyne phase-locked loop with 
a post-loop filter F{u) that realizes optimal smoothing. 

The variance of the optimal frequency-domain 
smoother is P, d, [T^ 



\Sxv{i^)['' 



Sy{uj) 



2fc(KA/fc-K 1)1/2' 



(4.22) 



which is the same as the smoothing error derived by the 
state- variable approach in Eq. (|4.5p . as expected. The 
interested reader is again referred to Ref. for an ex- 
cellent treatment of optimal frequency-domain filters and 
smoothers. 



V. QUANTUM POSITION AND MOMENTUM 
ESTIMATION BY OPTICAL PHASE 
MEASUREMENTS 

A. Quantum Kalman-Bucy filtering 

So far we have assumed that the mean phase of the 
optical field contains classical random processes to be es- 
timated in the presence of quantum optical noise. In 
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this section we investigate the estimation of inherently 
quantum processes carried by the optical phase. Specif- 
ically, we revisit the classic problem of quantum-limited 
mirror position and momentum estimation by optical 
phase measurements. First we review the problem of 
optimal real-time estimation by Kalman-Bucy filtering, 
previously studied by Belavkin and Staszewski [s^ and 
Doherty et al. J||. 



CW laser beam 



to PLL 




u 



FIG. 6: (Color online). Position and momentum estimation 
by optical phase measurements. 

We model the mirror as a harmonic oscillator, as de- 
picted in Fig. m 



dx p 

dt m 
dp 

'di 



muj.„x + 



2Mhu!o cos 9 



(5.1) 
(5.2) 



where x{t) and p{t) are quantum position and momen- 
tum operators, m is the harmonic-oscillator mass, 
is the mechanical harmonic-oscillator frequency, the last 
term of Eq. (|5.2p is the radiation pressure term, M is the 
number of times the optical beam hits the mirror, 9 is 
the angle at which the optical beam hits the mirror, and 
I{t) is the optical flux operator, consisting of a mean and 
a quantum noise term. 



V 

m = T- 

nujo 



Aiit). 



(5.3) 



The constant force term can be eliminated by redefining 
the position of the harmonic oscillator, so we shall neglect 
the constant radiation pressure term from now on. AI{t) 
is approximately a white Gaussian noise term for a high- 
power optical coherent state. 



(5.4) 



A/(t)A/(r) 5(t- 



The Gaussian approximation neglects the discreteness of 
photon number, and is valid if the number of photons 
within the relaxation time of the filter impulse response 
is much larger than 1. We can then write the quantum 
system model as 



d 


X 







1/m' 




X 




■q" 


dt 


P 











p 


+ 


1 



(5.5) 



with initial conditions 

(Hto)) = 0, 



mo)) = 0, 



(5.6) 



and radiation pressure acting as the quantum Langevin 
noise. 



U = 



{u{t)u{T)) ^US{t-~T), (5.7) 
-- — , (3=2Mkocos9. (5.8) 



The mirror position is observed via optical phase mea- 
surements using a phase-locked loop. In the linearized 
regime, we can define the quantum observation process 
as 



y^[P 0] 

(z{t)z{T)) = ZS{t- t), Z : 



(5.9) 
(5.10) 



Our linearized model is consistent with the general model 
of continuous quantum non-demolition (QND) measure- 
ments (8i, i9i, il6]. 

To apply Kalman-Bucy filtering to the estimation of 
mirror position and momentum, let us define 



Aa; = X — x\ Ap = p — p , 
(A£2) I 



(5.11) 



{AxAp + ApAx) 



'AxAp + ApAx) 
(Ap2) 



(5.12) 



In the linearized model, the Wigner distribution re- 
mains Gaussian and non-negative provided that the ini- 
tial Wigner distribution is Gaussian, so it can be regarded 
as a classical phase-space probability distribution, x[t), 
p{t), AI{t), and y{t) can be regarded as classical ran- 
dom processes with statistics governed by the Wigner 
distribution, and we can apply classical estimation the- 
ory directly. The off-diagonal components of the variance 
matrix are written in terms of symmetrized operators to 
ensure that they are Hermitian and also obey Wigner- 
distribution statistics. 

The Kalman-Bucy variance equations hence become 



dt 
dt 

dT.22 

dt 

V : 



1 

TO 
1 



(S]2i + Si2)-FI]?i, (5.13) 

-E22 'Tiw^Eii - ySiii;i2, (5.14) 



-TOCj2„(Si2 + E21) - ySi2S21 -f f/, (5.15) 

4/327? 



hujQ 



(5.16) 



The steady state is given by the condition d'S/dt — 0. 
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After some algebra, 
h ^/2 



(511) s 

(512) s 
(S22)s 

where 



2mu!m Q 

(S2l)ss — 



h 1 

2Q 



1 + 

1 - 



Q 



(i + Q 



2^l/2 _ ^ 

2\l/2 , 



1/2 

- 1 
1/2 



1 + 



Q = 



'UV 



(5.17) 
(5.18) 

2)1/2^ 

(5.19) 
(5.20) 



is a dimensionless parameter that characterizes the 
strength of the measurements. The position uncertainty 
(Ai^)ss = (Sii)ss is squeezed due to the continuous 
QND measurements, while the momentum uncertainty 
(Ap^)ss = (S22)ss is anti-squeezed due to the radiation 
pressure. These results have also been derived by various 
groups of people [§, '9^1 , although here we have shown how 
one can realistically implement the optical measurements 
of a mechanical oscillator. 
The Kalman-Bucy gain is 



Apr 



Sii 
S21 



(5.21) 



The estimator equation becomes 



d 


' x'' 




-VT,ii l/m' 




' x'' 


di 


.P' . 




-yE2i-mw2„ 




P' 



V 



Sii 
S21 



y, 



(5.22) 



where y is the measurement record of y. The filter relax- 
ation time is on the order of 



1 



1 



nSii)ss V2a;„[(l + 02)i/2_i]i/2' 



(5.23) 



which decreases for increasing Q, so the steady state can 
be reached faster for a larger Q. For Q ^ 0, — > 00, 
and a steady state does not exist. The photon number 
within the filter relaxation time is much larger than 1, 
and the assumption of white Gaussian radiation pressure 
noise is valid, when 



Q 



Vtf 

' 2V2/32n[(i + 02)1/2-1]! 



/2 



> 1. 



(5.24) 



On the other hand, the threshold constraint, which en- 
sures that the linearized analysis of the phase-locked loop 
is valid, is 



/3'(Si 



[3-'h [(i + o')'/'-i] 



1/2 



V2; 



miOri 



Q 



<1. (5.25) 



This condition, apart from a factor of 4, is the same as 
the large-photon- number assumption given by Eq. (|5.24p . 
and ensures that the linearized system and measurement 
model is self-consistent. 

The mirror position-momentum uncertainty product 
at steady state is 



{Ax%s{^f)ss = (S11S22) 

~ To^ 



1 + 0'- (1 + 



2\l/2' 



(5.26) 



and satisfies the Heisenberg uncertainty principle for all 
0, as one would expect. Furthermore, the covariances 
satisfy the following relation for pure Gaussian states [1, 



det(S,,) = (EnS22 - S?^) = -, (5.27) 



4 ' 



indicating that the harmonic oscillator conditioned upon 
the real-time measurement record becomes a pure Gaus- 
sian state at steady state. 



B. Smoothing errors 

From the classical estimation theory perspective, we 
should be able to improve upon Kalman-Bucy filtering if 
we allow delay in the estimation and apply smoothing. 
Here we calculate the steady-state smoothing errors us- 
ing the two-filter approach described in Sec. IIVBI The 
steady-state smoothing covariance matrix is 



1 I D-IN 



(5.28) 



where Sgs is the steady-state forward-filter covariance 
matrix, already solved and given by Eqs. (|5.17p - (|5.19p . 
The backward-filter covariances obey the following equa- 
tions: 



(5.29) 
(5.30) 





1 


dt 


m 


dEi2 


I 


dt 


m 


d'B.22 




dt 


— — r 



The steady-state values for the backward filter turn out 
to be almost identical to the ones for the forward filter, 

(■^ll)ss = (Sll)ss) (■^12)38 — — (Si2)ssj 

(S22)ss = (S22)ss, (5.32) 



and also satisfy the pure-Gaussian-state relation 

det(H..) = (SnS22 - S?^)^^ = ^. (5.33) 
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After some algebra, 



(nil),. 


h 




1 


1 








^ (l-^Q)l/2 


(ni2)ss 










(n22)ss 


hmujm 
8 









(5.34) 
(5.35) 



These results can be confirmed using the frequency- 
domain approach outlined in Sec. IIV CI The position- 
momentum uncertainty product becomes 



(nii)ss(n22)ss — ^ 



1 



(l+Q2)l/2 



(5.37) 



uncertainty principle would allow, quantum mechanics 
seem to forbid one from experimentally verifying the vi- 
olation. In this sense the apparent paradox is analogous 
to the Einstein-Podolsky-Rosen paradox [it'I and may 
yet have implications for the interpretation of quantum 
mechanics. 

In practice, while one may argue from a frequentist 
point of view that delayed estimation of quantum pro- 
cesses is meaningless if it cannot be verified, smoothing 
should still be able to improve the estimation of a clas- 
sical random process in a quantum system, such as a 
classical force acting on a quantum harmonic oscillator 

M. 



which is smaller than the Heisenberg uncertainty product 
/A by 4 to 8 times. 



VI. CONCLUSION 



C. Discussion 

While counter-intuitive, the sub- Heisenberg uncertain- 
ties given by Eqs. (|5.34p - (|5.37p do not violate any ba- 
sic law of quantum mechanics. The reason is that we 
only estimate the position and momentum of the mir- 
ror some time in the past as if they were classical ran- 
dom processes with Wigner-distribution statistics, but it 
is impossible to verify our estimates by comparing them 
against the mirror in the past, which has since been ir- 
reversibly perturbed by the unknown radiation pressure 
noise. In classical estimation, x{t) and p{t) are classical 
random processes unknown to the observer but can in 
principle be perfectly measured or simply decided at will 
by another party, so it is possible to compare one's de- 
layed estimates against the perfect versions and verify the 
smoothing errors. In the quantum regime, however, one 
cannot measure the mirror in the past more accurately 
without disturbing it further, and the only way for us 
to obtain perfect information about the mirror position 
or momentum is to perform a strong projective measure- 
ment. Unlike the Kalman-Bucy estimates, which predict 
the mirror position and momentum at present and can 
be verified by performing a projective measurement at 
present, it is obviously impossible to go back to the past 
and perform a strong projective measurement on the mir- 
ror to verify our delayed estimates without changing our 
model of the problem. 

It is also impossible to perfectly reverse the dynamics 
of the mirror in time and recreate the past quantum state 
without introducing additional noise, because quantum- 
limited optical phase measurements prevent us from ob- 
taining any information about the optical power fiuctu- 
ations, and the dynamics of the mirror subject to the 
unknown radiation pressure noise is irreversible. Thus, 
even though classical estimation theory indicates that we 
can achieve more accurate estimates than the Heisenberg 



In conclusion, we have used classical estimation the- 
ory to design homodyne phase-locked loop for quantum 
optical phase estimation, and shown that the estimation 
performance can be improved when delay is permitted 
and smoothing is applied. We have focused on coher- 
ent states, as it can be regarded as a classical field with 
additive phase-insensitive noise upon homodyne detec- 
tion, and classical estimation techniques can be applied 
directly. The optimal adaptive homodyne measurement 
scheme for nonclassical states remains an open problem. 
Along this direction Berry and Wiseman have recently 
suggested the use of Bayesian estimation for narrowband 
squeezed states when the mean phase is a Wiener pro- 
cess fT§|. Generalization of their scheme to more general 
random processes is challenging but may be facilitated 
by classical nonlinear estimation techniques [l, 0| ■ 

When we apply the same classical techniques to the 
quantum-limited estimation of harmonic-oscillator posi- 
tion and momentum, we find that the two conjugate vari- 
ables can be simultaneously estimated with inferred ac- 
curacies beyond the Heisenberg uncertainty relation, if 
smoothing is performed. Although quantum mechanics 
seems to forbid one from verifying the delayed estimates 
by destroying the evidence, this result remains counter- 
intuitive and may have implications for the interpretation 
of quantum mechanics. In the general context of quan- 
tum trajectory theory [23|, quantum smoothing deserves 
further investigation and should be useful for quantum 
sensing and communication applications [H, [2l| . 
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